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ABSTRACT 


Improvement B to the partial-ref lection-drifts experiment at Urbana, 
Illinois are completed. The results of the improvements include real-time 
processing and simultaneous measurements of the D region with coherent 
scatter. Preliminary results indicate a positive correlation between drift 
velocities calculated by both methods during a twc-day interval. 

The possibility now exists for extended observations between partial 
reflection and coherent scatter. In addition, preliminary measurements 
could be performed between partial reflection and meteor radar to complete 
a comparison of methods used to deteruiine velocities in the D region. 
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1 . INTRODUCTION AND STATEMENT OF THE PROBLEM 

The purpose of this thesis is to describe how the D region of the 
atmosphere is remotely probed by a ground-based radar system used at the 
Aeronomy Field Station near Urbana. The method used is referred to as the 
partial-reflection-drifts experiment. Other experiments which use different 
methods at this installation include coherent scatter and meteor radar. 

The partial-reflection-drifts experiment is used to determine the 
horizontal-drift velocity of ionized irregularities in the ionosphere. If a 
point radio source is used, the horizontally stratified irregularities 
produce a diffraction pattern over the ground. By sensing this diffraction 
pattern with a minimum of three antennas the horizontal-drift velocity can be 
computed. 

The former system at the Aeronomy Field Station lacks the ability to 
produce reasonable amounts of data so that the dynamics of the upper atmo- 
sphere can be studied. The amount of computer time needed to process data 
is one of the limiting factors of the former system. Since the computer 
used for the drifts experiment also is used by the other experiments, 
simultaneous measurements for comparing techniques are difficult. The system 
also lacks a suitable storage medium for the results so that advanced studies 
could be performed. 

The new system solves the problem by modifying the way data are col- 
lected, processed, and stored. By using small inexpensive microcomputers to 
distribute the computation load, simultaneous measurements become possible. 
The postprocessing computation time can be reduced by performing some 
calculations during the collection of data. The ability to store results on 
a floppy disk makes the storage medium ideal for postanalysis of results. 
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The mejor objective of these enhancementt it to optinice the usefulness 
of the partial-ref lection-drifts technique as a valid method for studying 
the dynamics of the upper atmosphere. 
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2. THEORY OF PARTIAL-REFLECTION DRIFTS 

To determine the hprl^fJhtal-drift velocity of the ionotphAre it ip 
neceiiary to illuminate the ionoaphere vith a aingle radio-wave point source. 
When this ia done a diffraction pattern is formed from the ionised irregu- 
larities in the D region. According to Erigga (1977) the ions in the region 
of interest move vith the neutral air« the drift velocity of the ions can be 
measured by applying the method developed by Briggs et al.. (1950) known as 
the full-correlation analysis. The radio waves are directed at vertical 
incidence, it is therefore necessary to have a system of spaced antennas to 
track the motion of the irregularities. With three spaced antennas to sample 
the amplitude of the diffraction pattern observed on the ground, spatial 
properties and the movement of the pattern can be deduced. 

The diffraction pattern is sampled at all antennas at equally spaced 
time intervals. The resulting records of amplitude represent the fluctua- 
tion of the diffraction pattern as it moves across the ground. The method of 
analysis applied to the recorded amplitudes requires the autocorrelation for 
each observation point and cross correlation for every pair of observation 
points. It is assumed that the information contained in these correlation 
functions can completely determine the velocity, direction of drift, the siae 
and the orientation of the irregularity which is inferred from contours of 
constant correlation which are in the form of ellipses (Briggs 1968). 

The signal amplitude of the diffraction pattern can be represented by a 
function R(x,y,k), where x and y represent the space coordinates and k the 
time coordinate. With the amplitudes of the diffraction pattern defined by 
a function of three variables, two of them being spatial variables, the 
corresponding correlation function is also a function of space and time. 
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For the practical analyaia the apatial variablea are renoved and the 
diatancea and direction of the baaelinea between antenna* are recorded for 
the latter part* of the analyaia. By numbering the antenna ayatem aa ahown 
in Figure /?.!• the diacrete correlation function* are function* of tine only 
and have the general forma; 
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N-T 
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for the autocorrelatin&c» and 
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for the cross correlations. 

The computed autocorrelations from each antenna should be identical, 
but due to statistical variations they are not in practice and the mean auto- 
correlation is used for the analysis. The cross correlation function has 
a maximum value displaced from the origin by some amount t' . The greater 
the velocity of the irregularities the smaller t* is. It is important, 
therefore, to separate the antennas by a sufficient distance so that t' is 
measurable, but not to space them so far apart that the cross correlation is 
zero (Briggs 1976). For the frequency used for our antenna system the optimum 
separation for the antennas is about 160 meters for D-region measurements. 

The antenna separation for the Urbana drifts experiment is 169 meters for the 
shorter sides of the triangle and 240 meters for the hypotenuse side. 

The typical shape of the resulting auto- and cross correlation func- 
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Figure 2.1 Arrangement of receiving antennas. 
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tiont reiemble* thoae in Figure 2.2. 

The inforaation that is extracted from the correlation £unctioi&, for the 
full correlation analysis is; 

1) The maximum value of the cross correlation (p ), , and the asso- 

m ij 

dated time displacement 

2) The values of time displacement autocorrelation 

that correspond to the maximum values (p ). . found for the cross correlation. 

m ij 

The method that was first used for analyzing drifts records was known 

as the method of similar fades (Mitra 1949). This was a simple approach to 

the problem of determining the drift velocity. The major disadvantage of 

using this method is that it assumed the irregularity did not change in 

shape as it moved across the observation points. The full correlation 

analysis corrects this assumption by computing the correlation contour that 

represents the fading of the drifting pattern with no random changes. It is 

assumed that this contour is in the shape of an ellipse but there is no 

physical reason why this should be so, and it is possible that the contour 

might be of any shape whatever (Briggs et al., 1950). To calculate the 

ellipse the values of (t )^ , and (t*),, are used in the formula; 

m ij ij 

(t )^ = (t') ^ - (t ) ^ 

where corresponds to the time shift for the autocorrelation to tal- 

on the value of the cross correlation at zero lag. By using the formula 
above the values of t^^j are less affected by statistical variation because 
it is computed at a higher level on the correlation functions (Briggs 1957). 
Three velocity vectors defined by 

magIv^ 

' ■'.J ' 3 

0 . 


'Ij 

ANGLE I V 




ij 


ij 
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Figure 2.2 


(a) Autocorrelation function . 


(b) Cross correlation function. 
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are computed. The ' a are the baaeline distances in meters and the 
are the angles measured from receiver "i" to receiver with respect to 
the north. The ellipse that defines this contour of constant correlation is 
referred to as the characteristic ellipse. 

The characteristic ellipse has the dimensions of velocity. If the 
ellipse is multipled by the half time of the autocorrelation, the resulting 
ellipse will have dimensions of meters and gives the relative size of the 
irregularity. The ellipse is rotated through some angle 6^ and has semimajor 
axis a and semiminor axis h. The axial ratio a/b gives an indication of the 
elongation of the irregularity and the angle 6^ gives the direction of 
elongation. 

The next calculation necessary is to determine the apparent velocity 
components between the antenna pairs. The three vectors 


V' ■ 

12 


®12 

V “ 

13 

d / T ' 
13' ^13 

®13 

V ■ 

23 

d /t* 
23' 23 

®23 


are computed and the end points will lie on a straight line referred to as 
the V'-line. If a vector perpendicular to this line is drawn from the 
origin it will define the magnitude and direction of the apparent velocity. 
The important information that is used from this is the slope of the V'-line 
defined by these vectors. The point of tangency of this line to the ellipse 
gives the true direction of drift <P measured clockwise from north. The 
vector from the origin to the tangency point is known as the vector. 


The true velocity is then given by; 

(V ) ^ 

V 


drift V 

This velocity must be divided by a factor of two to yield the correct 
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result* due to the point source effect. The east and north components of 
the drift are given by; 

''e ■ ''n ■ 

The geometry of the full-correlation analysis is shown in Figure 2.3. 

The full-correlation analysis has been used for several years. During 
that time many investigators in the field have done extensive studies to 
determine the validity of the analysisi by producing models of the 
ionosphere and performing computer siiitulations (Fitteway et al., 1971). The 
best way to determine the validity of this analysis is to compare results 
from simultaneous measurements. This is done in Chapter 5. 
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3. SYSTEM DESCRIPTION 

3.1 FORMER SYSTEM 

The partial-reflection system supports two experiments that are useful 
for D-region studies. The experiment that is being enhanced is the drifts 
experiment, the other is the differential-absorption experiment which is 
used to measure electron-density concentrations. A block diagram of the 
former Urbana partial-reflection system is shown in Figure 3.1. The major 
elements of this sytem are the transmitter, transmitting and receiving 
antennas, receiver, data-acquisition system, and radar controller. 

The transmitter is a multistage tube type and is fully described by 
Henry (1966) and by Pirnat and Bowhill (1968). The final output is split 
between two 50 ohm coaxial cables. This is done to insert a phase shift of 
90 degrees so that when the two outputs are fed into an orthogonal-dipole 
array the transmitted wave will be circularly polarized. Attenuators are 
used after the phase shift to ensure circular polarization. The characters 
of this transmitter are: 

Peak Power 35 kW 

Frequency 2.66 MHz 

Pulse Width 23 psec 

Output Impedance 50 ohm, unbalanced 

Pulse Repetition Rate 10 Hz 

The layout of the antenna system is shown in Figure 3.2 (Weiland and 
Bowhill 1981). The array west of the Field Station Building is the trans- 
mitting array and the array to the east is the receiving array. Both antennas 
were identical when they were first constructed. The receiving array has 
been divided into four quadrants so that the drifts experiment can be 
implemented. 

The transmitting array consists of 60 half-wave dipoles. The dipoles 
in the N-S and E-W direction are coupled together to form an orthogonal set. 
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Figure 3.1 Block diagram of the former partial-reflection system. 
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This i* done to achieve circulatly-polarixed-tranamitted waves. By inserting 
a relay network that inserts a 180 deg phase shift, the ordinary or extra- 
ordinary pulse can be transmitted. For the drifts experiment only one 
polarization is necessary and the ordinary is perferred over the extraordi- 


nary because the ordinary mode is attenuated less in the D region. 

The receiving array consists of four smaller arrays called quadrants, 
that contain three full wavelength dipoles in both the N-S and E-W direction. 
The N-S and E-W dipoles are separately matched to one 30 ohm coaxial cable 
that is fed into the Field Station Building. The selection of any quadrant 


is accomplished by a relay network mounted on the inside east wall of the 
Field Station Building. 

The receiver used for partial reflection is extremely linear. This is 
a requirement because the signal amplitudes must be accurately measured over 


a 50 dB range (Weiland and Bowhill 1978). The characteristics of the recei- 


vers are listed below: 


Center Frequency 
Noise Figure 
Bandwidth 
Recovery Time 

Gain Variation 
Rf Input Impedance 
Output Impedance 
Linearity 

The data-acquisition system for 
Digital Equipment Corporation PDF 15 


2.66 MHz 
3 dB 
40 kHz 

200 ysec after removal of .IV 
RMS input 
3 dB max 

50 ohm, unbalanced 
10,000 ohms 

55 dB for 1 dB deviation 
former system consisted of a 
computer for processing of data and 


a Hewlett-Packard 5610A analog-to-digital converter for digitizing data from 


the receiver. The output of the A/D was a 10-bit word in two's complement 
form. The entire range of the A/D was not used however since the input 
range for the A/D was 1 volt and the output of the receiver is zero to -i-1 
volt. The conversion rate is 100 kHz and corresponds to a height resolution 
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of 1.5 km. 

The radar controller referred to as the pulser generates the necessary 
control signals to operate the system. The pulser generated a 25 ’-micro second 
pulse to turn on the transmitter* a 100-microsecond long pulse to blank the 
receiver, a trigger pulae to activate the A/D and a logic signal to switch 
between the ordinary and extraordinary modes. 

3.2 SYSTEM MODIFICATIONS 

The purpose of modifying the system yas to improve its performance 
and adapt the sytem to work with a microcomputer. The goals of these 
improvements are to: 

1) Increase amount of data that can be taken daily. 

2) Decrease postprocessing time. 

3) Enable simultaneous measurements between the drifts experiment 
and the coherent-scatter or meteor-radar experiment. 

4) Form a data base of results for advance analysis. 

A block diagram of the new system is shown in Figure 3.3. The compon- 
ents of the former system that were replaced include the computer, pulser, 
and the A/D converter. The new system does not require both polarizations for 
collecting drift data. The new system therefore does not control the phase 
shifting networks for the extraordinary mode. The ordinary mode is prefer- 
red. The new pulser increases the pulse repetition frequency (FRF) from 
10 Hz to 200 Hz. It was considered that by increasing the FRF to this rate 
the signal could be oversampled and would improve the signs 1-to-noise ratio. 
The Apple II computer is being used to gain independence from other experi- 
ments that use the FDF-15 computer. A new A/D converter is also necessary 
for compatibility with the Apple computer. With partial-reflection drifts 
independent from other experiments, simultaneous measurements can be per- 
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Figure 3.3 Block diagram of new partial-reflection system. 
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formed. The organlBAtion of the computer ayitem will accomplieh the rest of 
the desired goals. 

3.2.1 Computer system . A block diagram for this computer system is 
shown in Figure 3»4. The computer system used for the drifts experiment 
consists of an Apple II Plus for the collection of data, an Apple II Plus for 
processing of data, a 10 megabyte Winchester disk drive for data storage, and 
various peripherals for the Apple XI computer. The peripherals consist of a 
real time clock, a parallel interface, a modem, and an arithmetic processor. 
The function of each of the peripherals is summarized below. 

1) Real-Time Clock. This peripheral generates the necessary informa- 
tion to record the date and time with each record of collected data. 

2) Parallel Interface. This peripheral forms the interface to the A/D 
converter and the relay network for quadrant switching. It contains four 
parallel ports that can be used for input or output. It also has the 
ability to generate square wave signals that are used to synchronize the 
collection program with the A/D converter. The main component on this 
interface card is the 6522 versatile-interf ace adapter. It is an all pur- 
pose parallel-interface adapter for the 6502 eight bit y processor chip 
which is the Apple's CPU. 

3) Arithmetic Processor. This powerful peripheral enhances the com- 
putation power of the Apple. It is used to perform the necessary multipli- 
cations for the correlation functions. The processor card uses the AMD 9511 
arithmetic processor chip. This integrated circuit can perform a variety of 
operations including transcendental functions. The time required to perform 
a 32-bit by 32-bit multiplication is 105 microseconds. 

4) Modem, This peripheral along with the necessary software enables 
data transfers between computers. This modem is a direct-connect type with 
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Figure 3.4 Block diagram of new computer system. 
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a tranafer rate of 300 baud. 

3.2.2 A/D interface . Thia interface performs the necessary analog- 
to-digital conversion for the drifts experiment. The interface converts the 
amplitude detected signal from the receiver to an eight~bit unsigned binary 
number. This binary number is transferred into the computer by means of a 
parallel interface. The circuit diagram for the interface is shown in 
Figure 3.5. 

The A/D performs 10 ^ conversions per second. This corresponds to a 
range resolution of 1.5 km. The main component of the interface is XC8. 

This integrated circuit package is a successive approximation type A/D with 
a conversion speed of six microseconds. For a sampling rate of 100 kHz 
that leaves four microseconds to sample the signal and latch the data to 
the parallel interface. 

Only one control signal is needed to drive the interface and is gene- 
rated by the parallel interface by means of an internal timer. This signal 

is a symmetric-square wave with a frequency of 100 kHz. The purprse of 1C4 

is to produce an asymmetric-square wave since the sampling and conversion 
time cannot be equal for the A/D used. 

The operational amplifiers ICl and IC2 are used to buffer the input, 
amplify, and level shift of the signal if necessary to accommodate the analog- 
to-digital converter. For the A/D used the level shifting ability is not used 
but the signal must be amplified by a factor of ten to take advantage of the 
range of the A/D. Because damage may occur if the input to 1C8 is greater 

than ■i'lO volts or less than zero, Dl and DZ are used to limit the signal to 

this range. The free running voltage-controlled oscillator (VCO) chip IC6 
is used as the clock for the A/D and is buffered by 1C7. The clock frequency 
for a six-microsecond conversion time is 1.33 MHz. The digital outputs of 
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the A/D ere buffered by 1C9, ICIO and 1C5. 

3.2.3 Pulaer . The pulaer the main controlling unit of the drifts 
experiment. This unit produces the necessary control signals for the trans- 
mitter! receiver blanker, and computer. The timing diagram that shows the 
relationship between these signals is shown in Figure 3.6, and the circuit 
diagram is shown in Figure 3.7. 

All the control signals of the pulser are derived from ICl which is a 
VCO. The resulting signal from this chip is fed through divide-by-two and 
divide-by-ten counters to produce five separate FRF rates. The desired FRF 
rate is used to drive four monostable circuits which produce the control 
signals. 

The ability to choose the FRF rate allows some flexibility for future 
work with the drifts or the electron density experiment. The pulse width of 
the transmitted pulse can easily be varied by adjusting the potentiometer on 
IC4. When adjusting the FRF and pulse width the maximum duty cycle of the 
final stage amplifier of the transmitter must not be exceeded. 

The blanker signal is set high before, during, and after the trans- 
mitter pulse is activated. The interrupt request control signal for the 
computer is adjustable by means of a panel mount potentiometer. The adjust- 
ment of this signal determines the approximate starting height where data are 
collected by the computer. 

Due to the problems observed when both transmitters were operating for 
partial reflection and coherent scatter, the pulser was modified. The FRF of 
the coherent- scatter transmitter is 400 Hz and the TTL level signal that was 
generated by the radar controller for coherent scatter was used in place of 
ICl. The circuits that generate various FRF rates were disconnected until 
the need for them arises. 
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Figure 3.6 Timing diagram for new pulser. 
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4. DATA COLLECTION AMD ANALYSIS SYSTEM 

4.1 OBJECTIVES OF THE ANALYSIS 

The major objective of the analysis is to efficiently compute the drift 
velocity. The analysis can be broken down into two phases. The first phase 
is the collection of data and the second is the postprocessing of those data. 

Since the actual collection of data requires very little time, it is 
advantageous to use the idle time of the computer to start processing data as 
they become available. This eliminates two problems in the procedure. The 
first is the amount of storage required to store the actual amplitudes, and 
the second is the time required for postprocessing. By computing the 
correlation functions at the time of the collection, the amount of memory 
required is reduced by a factor of two. The amount of postprocessing time 
is dramatically reduced since multiplication is a very time consuming 
operation. 

The first phase now accomplishes more while the second phase is 
required to do less. This will optimize the procedure of the analysis in the 
sense that all usable time in the collection phase is used to reduce memory 
constraints and postprocessing time. 

4.2 COLLECTION SOFTWARE 

The collection software for the drifts experiment is written primarily 
in machine code. This is done to optimize the performance of the program so 
that it will execute at the highest speed possible on the computer. The 
outer interpreter language used to link the machine code to the computers 
operating system is FORTH. This higher-level language is extremely flexible 
and fast. 

The collection program performs two major functions. One function is to 
collect data from the A/D and control the quadrant switching network. The 
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other is to compute the auto- and cross correlations. During each transmit 
pulse the radar controller interrupts the processor to collect data. The 
routine that is executed when this occurs is called the interrupt service 
routine (ISR). While the processor is waiting for the next interrupt signal 
from the radar controller it computes the correlation functions. This is the 
second function of the software. 

The flow charts for the ISR are shown in Figures 4.1a, 4.1b, and 4.1c. 
This section of the program is executed after every transmit pulse. The 
routine can be in one of three states during execution. 

The first state is called the stall state (Figure 4.1c). Because 
switching between quadrants is done with relays, sufficient time must be given 
to allow the relays to switch. This segment of the program counts s specified 
number of interrupts for this delay, k stall counter and a stall flag control 
the flow of this state. If the flag is set to "1" the program counts the 
interrupts until the maximum count is reached. For the present version of 
the collection program the maximum count is ten. This gives the relays 50 
millisecond to switch for a FRF of 200 Hz. 

The second state is called the wait state (Figure 4.1b). During this 
time the stall flag is set and the next antenna is switched into the receiver. 
One of the output ports on the parallel interface controls the quad select. 

By writing the appropriate number to the J2 port one of the four antennas is 
sampled; the following list shows the number and the quadrant it selects. 

QUADRANT VALUE IN HEX TO SELECT 

NE OE 

NW OD 

SW OB 

SE 07 

After all quadrants are sampled the pointers in the program are updated. 
The sample ready flag is set for the processing part of the program. The 
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Figure 4.1a Flowchart for interrupt service routine (ISR) 
collection state. 
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Figure 4.1b Flowchart for (ISR) wait state 
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values for the nev sample are summed to calculate the mea^\ signal level for 
each antenna at each altitude. If enough samples have been collected for 
one minute of data, the end-of-record flag is set, and all other flags arc 
cleared. If not, the program checks to see if the stall count is set for "1" 
and either increments the stall counter or clears the stall flag. 

The last state of the ISR collects the data from the A/D (Figure 4.1a). 
The first step during this phase is to synchronise the A/D start convert 
signal with the program. Timer 1 is used on the versatile interface adapter 
(VIA) for this control signal. It generates a square wave with a frequency 
equal to 10 /(2N 4) where N is the value in the timer 1 latch. The output 

port will toggle PB7, which is the MSB of the port, each time the timer 1 
counts down the value in the latch. By writing to the output port, FB7 is 
reset and the countdown restarts. Once the A/D is synchronised a dummy read 
is done to enable latching on the input port. 

The program also activates one of the annunciators on the Apple game 
l/O port by writing and reading a memory location. This gives an indication 
when the program is reading data. This signal which is generated on pin #15 
of the game port is used to set the altitude range gate. The program then 
reads 20 samples into memory and sums eight of those samples into the input 
buffer. This summation is performed ten times for each antenna and altitude. 
This oversampling improves the quantization of the signal and is recoismended 
when the A/D has a small word length (Hagen et al., 1973). 

When the computer is not busy with the ISR it is computing the correla- 
tion functions. The routine that performs this operation is called DCORR. 

The flowcharts for this routine are shown in Figure 4.2. The routine DCORR 
computes the auto- and cross correlations in this manner it is necesary to save 
the past "L" samples, where "L" corresponds to the number of lags that must be 



Figure 4.2 Flowchart for main code DCORR. 
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conputed. By noting that for a croaa correlation Rab (T) where A and B are 
the two different aeriea, forma the negative lags for R^(T). For the 

four-antenna syatem, four autocorrelationa and six croaa correlations must be 
computed. The collection program computes these by comput’ng sixteen cor- 
relations with five positive lags and the zero lagi This simplifies the 
program and reduces the execution time to compute the correlation functions. 

After each sample becomes available for processing the program adds the 
products that have been computed with the past L samples to the lag table. 
This operation can easily be understood if the products added to the lag 
table are displayed in the format shown in Figure 4.3. In this example the 
first two lags plus zero lag are computed. At a time Tl the first sample is 
available for processing. The products formed at this time interval are 
added to the corresponding lag bins. At the next time interval the calcula- 
tion is repeated with the lag reference advanced by one. By repeating this 
process the values in the lag bins will contain the desired results at the 
end of the sampling interval. For the drif ts-collection progran the process 
is repeated 128 times, this corresponds to a 51-second sampling interval. 

4.3 POSTPROCESSING SOFTWARE 

The postprocessing software performs the full correlation analysis. The 
input to this program is the correlation functions computed by the collection 
software. The output generated includes various parameters of the analysis 
that are of interest. 

The flowcharts for the postprocessing program are shown in Figures 4.4a, 
4.4b, 4.4c, and 4.4d. Two different versions of this program have been 
written. The first version, which follows the flowcharts, calculates the 
drifts velocity four times with different combinations of three antennas. 
Because of an unexpected failure in the relay network for the NE quadrant, a 
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Figure 4.3 Method used to calculate correlation functions (refer to text). 













Figure 4i4a (coivt.) 
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Figure 4.4a (cont.) 
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Figure A. 4b Flowchart for computation of apparent and true velocity. 
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second version of the processing progran was developed. The second version is 
similar to the first except that it is optimized to calculate the drift 
velocity once with one combination of three antennas. 

The major computations for the postprocessing program are: 

1) Complete the calculation of the auto- and cross correlations. 

2) Interpolate the cross correlation to datermine the true maximum. 

3) Compute the ellipse and its parameters. 

4) Compute the apparent and true velocities. 

The collection software does not have time to complete the calculation 
of the correlation functions. It stores all the necessary infoirmation so that 
the postprocessing program can complete the computation. The actual data 
that the processing program has to work with are the sum of the products and 
the mean values. The sums of the products^ which are more commonly referred 
to as the autocorrelation and cross correlation, are four-byte integer values. 
The autocovariance is computed by subtracting the square of the mean and the 
cross covariance is computed by subtracting the product of the means. After 
the cross covariances are normalized using the zero lag of the autocovariance 
the autocovariance is normalized. The use of the terms autocorrelation and 
cross correlation throughout the rest of the text refer to the normalized 
autocovariance and the cross covariance. 

Once the correlation functions are computed the next calculation is to 
interpolate the maximum of the cross correlations. For this interpolation a 
polynomial of the fourth degree is used. Using Stirling's Formula the first, 
second, third and fourth order differences are computed. If the maximum value 
is denoted by y^ the required differences are 

Ay_i = - y_i = yi - 3y^ + 3y_^ - y_2 

= y2 - 3y^ - y_^ 

A2y_i = y^ - 2y^ + y_^ = y2 " + ^y^ - 4y_j^ + 
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The form of the polynooiel is: 

y " •<, + a^x + X2 x2 ♦ «3x 3 + e^x^ 

The values of the "a" constants for this equation are: 

% - Vo 83 “ 1/12 (aV_2 

- l/2(Ay_^ + Ay^) ~ l/12(A3y_2 + A3y_^) - 1/24 aV_2 

82 ■ l/2(A2y_j^) - l/24A**y_2 

To determine where the maximum occurs it is necessary to set the deriva- 
tive equal to zero. 

^ + 2a2X H- Sa^x^ + ^8^x^ * 0 

The roots to this equation can be found by applying Newton's method. This 
requires the next derivative: 

d^V 2 

— r = a, + 2a„ + 6a,x + 12a, x 

dx'^ ^ 

The first approximation to the root is zero since the equation has been formed 

with respect to the maximum value y located at x > 0. By iteration the 

o 

interpolated value of x is found by using 

V = V - ^ 

^ Vl dx/ ^^2 

"N" times until the root converges to a specified tolerance. The value of N 
for the processing program is 100. Most of the time only a few iterations 
are necessary for convergence. 

Once the value of x is determined the maximum value is computed by 
evaluating the polynomial. The time delay in seconds is computed from 
TD - [(KM - 6) + x] *.4 + CA 

The value of KM is the index where the maximun occurred. For the zero lag the 
index KM would have the value of six. For positive lags take on the values 
7-11, for negative lags the values 5-1. The values of "CA" is an antenna 
cycling correction value since all antennas are not sampled at the same time. 



ORIGINAL PAGE IS 
OF POOR QUALITY 


44 


The approach taken to interpolating the nuiximum is a great inprovement 
over the former processing program. That program Was an adaption of the 
analysis progrsm used at the University of Adelaide (Australia) for drifts 
analysis (Welland in Edwards 1978). The former program did not have an 
iteration limit for the determination of the maximum, It also used a dif- 
ferent method for solving the maximum that wasted time, and at tiaies lacked 
accuracy. 

The second version of the processing program also did a parabolic inter- 
polation if the iteration limit was exceeded or if the maximum of the cross 
correlation was outside the range of plus or minus the third lag. Because the 
fourth order fit requires two values on either side of the maximum, and only 
the ±5 lags are Computed, the parabolic interpolation allows a maximum to be 
computed if it occurs on the ±4th lag. The fourth-order fit may not converge 
rapidly if the peak is broad. In this case a parabolic fit is the best 
approximation that can be made. The calculated value of x for this fit is 


l/2(- 


y-1 - y+i 


y_l - 2y^ + y + 1 

and the maximum correlation value is 


) 


V V 

■'max -^o 


(y-1 ~ y-n) ^ , (y-i 
2 


+ 2 
jj 


The next step in the calculation is the determination of the ellipse. 
From the analysis the ellipse is defined by the endpoints of the vectors; 
(Vc')j|^ 2 * general form of an ellipse in polar coordi- 

nates that has been rotated through some angle 6^ is 

coa^(® - 6o) , sln^(9 - 0o) 1 

2 ;:2 2 

a b r 

The three endpoints with polar coordinates 
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lufficient to deteraine the characterietics of the ellipae (Fooka 1965). 
By dsfining the tvo variablea: 




the value of Qq can be determined from the formula 


Dj^ 2 ('Cos 26 j^ - cos26^) - - cos26^) 

* Dj^3<sin20j^ - sin202) - - 810263) 


2 2 

The values of a and b are: 

2 (K3K2 - Kj_K^) 

a - 2 2 * 

(K2/r2^ - K^/r^ ) 

where 

1 + cos 2(01 - 60 ) 


V 1 + cos2(02 - Qq) 

*^3 ” 2 


2 (K 3 K 2 - Kj^K^) 

(Kj/r/ - K^/r/) 


V 1 - cos2(01 - 6 q) 
*^2 “ 2 

„ i - cos2(02 - 9 q) 
K-4 “ 2 


Because the correlation contours can have any shape whatsoever (Briggs 1950), 
the values of a^ and b^ may be negative. VHien this occurs the analysis is 
aborted. 

The last computation will calculate the apparent and true velocities. To 
calculate the apparent velocity the vectors, 73 ^ 2 ’ >^ 13 ' >^23* * computed. The 
magnitudes are 


V 

ij 43 

where d. . is the baseline distance between the antennas and t . ^ is the time 
ij ij 

displacement from the cross correlation. The angle for each vector is the 
angle measured clockwise from north from "i" to "j". The endpoints of these 
vectors define the V'-line. To compute the best fit, the method of weighted 
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least aquaret la used, the value o£ the W^’a is (r^) • By tranaforming the 
polar coordinatea o£ these endpoints to cartesian coordinatca (X,Y) the 
Y-intercept (YI) and slope of the line (M) are given by 

yi « EWYEWX^-EWXSWXY 
E W E WX^ - (E W X)^ 

EWEWXY-EWYEWX 

J»J m M .1— M l 

E W E WX^ - (E W X)^ 

Once the slope and intercept are calculated the apparent velocity can be com- 
puted. the apparent velocity is defined by the distance and direction of the 
line that joins the origin to the point (X,Y) perpendicular to the V'-line. 
the values of X and Y are: 


X = 


-m 


(m + _1 ) 
in 

and the apparent velocity in polar coordinates is 


(m^ + 1 ) 


V 




2.2 

X + y 


(j) » tan"^Cy/x) 

3 


The next step is to calculate the point of tangency between the ellipse and 
the V'-line. To calculate this point it is convenient to rotate the coordi- 
nate system so that cartesian coordinates can be used. The slope at any 
point (X^, Yq) on the ellipse becomes: 


m 




The slope m that is equated to this is 

-1 


m 




tan(<() - 0 ) 

^a o 


^ A 

2 W ^ 

a o 


The possible solutions to the tangency point are 

.2 


V = + 

X rv — - 


° (b^ + a^(m')^)^^^ 


X„ = ± 


2 , 

a m' 


° "(b^ + a^(m')^)^''^ 
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The correct valuei of end ere chosen based on the fact that the true 

velocity must be within 90 degrees of the apparent velocity (Fooks and Jones 

1961). The vector defined to the point (x^, y^) is known as the vector . 

* 

The vector that lies in the same direction as the (V^')y vector but extends 
to the V'-line is referred to as the V vector. To calculate this vector the 
y-intercept calculated for the V'-line must be rotated to the new coordinated 
system. The transformation is 

YI' ■ V^[sln((|i^ - e^) “ m' cos((|i^ - 0^)] 

If the slope of the V* vector is m^. The endpoint coordinates of the V' 
vector are: 


m^. - m' 



YI' 

- m' 


The magnitude of the true velocity is then 

V’ 

The direction of drift. 6^, must be modified by 6^ because the value calculated 
was in the rotated coordinate system. Therefore + 3^. 
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5. COMPARISOH OF PARTIAL-REFLECTION AND COHERENT-SCATTER WIUD DATA 

The new partial-reflection drifti experiment was used to make simul- 
taneous measurements with coherent scatter on April 24 and 25 of 1982 • These 
data are unique because during the time of the observations the upper 
atmosphere was exposed to energetic electrons from a solar flare. According 
to Briggs (1981), these two techniques used to determine atmospheric winds are 
basically the same, therefore the mean horizontal winds for the two days were 
compared between the experiments to determine if similarities existed. 

The horizontal-wind velocity computed by coherent scatter is a projection 
of the N-S and E-W components onto the line-of-sight direction for their 
antenna system. The antenna system is pointed in the south-east direction at 
an azimuth angle of 126 deg from the north, 1.5 deg off-vertical. The wind 
velocity components calculated by the partial-reflection-drifts experiment 
are projected onto the same azimuth for comparisons. 

Because Of the small angle from vertical that is used to determine the 
horizontal velocity for coherent scatter, hourly means are used for the 
comparison. The altitudes that overlap between the experiments are 70.5, 

75, 84, 88.5, 93, and 97.5 kilometers. 

The comparison of the results shows a positive correlation for five of 
the six altitudes and are shown in Figures 5.1 - 5.5. The points on the 
scatter plot should form a straight line with the slope equal to unity. The 
comparison at 79.5 km yields a correlation coefficient of .54 and a con- 
fidence level of 96.3%. The results at the next altitude had a correlation 
coefficient of .357 with a confidence level of 86.7%. These altitudes 
indicate a definite similarity between the two experiments. The comparison 
at 75 and 93 kilometers had a weak correlation. The correlation at 97.5 km 
was .524 with a confidence level of 95.8%. The slope of the scatter plot 
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Figure 5.1 Scatter plot at 75 km. 
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Figure 5.3 Scatter plot at 8A km. 
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Figure 5.5 Scatter plot at 97.5 km. 
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for this altitude does not match the expected slope for a one-to-one agree- 
ment. 

For some altitudes the agreement between coherent scatter and partial 
reflection is quite good. One observation about the altitudes that shows 
the greatest agreement is that they also correspond to regions where the 
scattered power is large, from 75 to 88.5 km altitude. The smaller correla- 
tions at 84 km and 75 km are somewhat unexpected if large scattered power 
improves the comparison. 

The computer simulation of the drifts experiments done by Pitteway et 
al. (1971) shows an interesting result. In the simulation the calculated 
drift velocity is plotted versus noise factor (Figure 5.6). The noise factor 
is a velocity perturbation, introduced to account for random changes in the 
ionosphere. If the calculated drift velocity tends to be larger than the 
actual drift, with an increasing noise factor, the scatter plots at 84 km and 
75 km may indicate this. By reducing the magnitude of the velocity calculated 
by the drifts experiment, better correlations result. The disagreement that 
occurs at 97.5 km may be due to receiver saturation. At this altitude the 
receiver saturates a small percentage of the time. During the collection of 
these data the objective of keeping the gain set for a range that is accept- 
able for all altitudes is hard to accomplish because of the large variations 
in the signal strength. The preliminary results from this simultaneous 
experiment indicate the need for future comparisons. 
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6. CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 

Tlie work has dealt with the enhancenents of the partial-reflection 
drifts experiment and the preliminary results of aimultaneoua measuremnnta. 
The modifications to the system: 

1) Allow simultaneous measurements with coherent scatter; 

2) Enable the drifts experiment to operate continuously so that large 
amounts of data could be collected; 

3) Decrease the time required for postprocessing of data which has 

a direct relationship to the amount of data that can be collected daily; and 

4) Permit easy access to processed results by storing them on floppy 

disk. 

The results of the preliminary simultaneous measurements have indicated 
an agreement between the drifts method and coherent-scatter's method for 
determining horizontal wind velocity. This set of results indicates the 
need to continue simultaneous measurements under careful control of the 
receiver gain, until enough data have been collected to assure statistical 
significance. 

Future improvements to the system could solve the remaining problems. 

1) ‘Signal-to-Noise-Ratio. By adding phase detectors so that coherent 
detection can be used the SNR would improve since the signal then can be 
integrated. 

2) Transmitter. If the transmitter power could be increased the lower 
limit of the altitude range could be reduced. 

3) Quadrant-Switching Network. This relay network should be replaced 
with electronic switching. The present design is not reliable for extended 


use. 
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4) Computer. If coherent detection is used, greater computing power 
is necessary to handle the complex correlation functions. 

5) Receiver Gain. The large variance in signal amplitude dictates 
that a programmable attenuator is necessary to attenuate signals at selected 
altitude range gates. 
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( RETURN PROM INTERN) PT ) 


THEN ROTATE POINTERS ) 

si!ST^I^8BSKSIs^''lif 1 


wpaw'iraPTML'M'iff 


( fix POINTERS POR ZERO LAG ) 


( SWITCH RELAY TO QUAOl FOR INPUT ) 


( SUN UP VALUES FOR THIS SAMPLE TO COMPUTE ) 
( MEAN OP EACH ANTENNAAND EACH HT ) 


i SeOO fy AD§' 

i s,'" 

OEY, 

if LDA, 

6 AD staJ 

1 AB sta; 

2 IF ♦ LDyJ 

3 So * LDA, 

4 BEGIN. 

I AD )Y STA, 

I AB )Y STA, 

7 DEY, 

LDA, 


( CLEAR NEXT LAC AREA FOR INCOMING VALUES ) 


SET FLAG FOR NEXT SAMPLE READY ) 


BD CMP, ( CHECK TO SEE IF ALL SAMPLES COLL 

1"^' INC, ! W Vc M"cS8WJ®B ?«B5 k 
OA • LDA, 


COLLECTED 


FOR F;U OF STALL PERIOD ) 


CO 

IFNE, 


J S Si” 

l?l 8S”( 

478 E4 

479 CO 


( RETURN PROM INTERRUPT ) 


( IF STALL COUNT - MAX COUNT CLEAR FLAGS ) 
( RETURN FROM INTERRUPT ) 


4 89 51 
490 BE 


( IF ALL SAMPLES TAKEN SET FLAG FOR END OF MINUTE ) 
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J701 

THEN, 

THEN, 


PL A, 

P; 

tax, 

LOA, 


( RETURN PRON INTERRUPT ) 


)DB 

j.L ISR 


R'EL A V ^ TO™ H A^ E , Q 
C^E sTj‘"'"’“ 




( END OF CODE FOR WAITISR ) 

FdIoSS?! tISeJI Bmwmi^AgfW 

CODE FOR STALL ISR } 


COONrSR''^POR StALL COUNT 


PLA, 

tat; 

PLA 


5^01 

ELSE. 


ik’, 

RTi; 


fOl 


1 ; 

1 : 
TAX, 
LDA, 
LDA, 
RTI, 

HIDc6pE 

:ISRCR^g^^ 

01A9 I C4C5 
LH , AAA5 , 
SC C, WAIT , 
80A9 , 8D c; 
C780 , 8D C, 
C058 , 8--' 
99 C,'4EI 
4EA5 4E^ 

8Eb1 

l8 BiB 
AD85 , AB 
HGATE_, 

? ?o8^* A 

18 C, ADB 

nn ; 

LOOP 
AAE6 


COMPARE FOR MAXIMUM COUNT ) 

IF <> THEN RETURN FROM INTERRUPT ) 


( OTHERWISE CLEAR COUNT AND STALL FLAG ) 
( RETURN FROM INTERRUPT ) 


( END OF CODE STALLISR ) 
INTERnjPT SERVICE ROUTINE 


48 C 
. 03D0 

,teo', 


ISR" ) 

6C C, ( CHECK FOR STALL AND INTEGRATION COUNT FLAGS ) 


63D0 



( RESET PB7 ) 

( READ IN DATA ) 


( ADD DATA TO BUFFER ) 


C701 


( RESTORE REG AND RETURN FROM ISR ) 


E8ES 

, A868 , AA68 

45A5 40 C, i 

: ISR J -5 DTOP +J 
ISRCR PLUSH DESTROY 
: PADDRC ^ WAITISR WAIT , . 

' STALLISR LW 1 ^ DCORR Cl + PADDR 
' DCORR HE + PADDR2 I ; 

i,8fM5o"?5 0 Kfflia‘>8!Lt“lrffi'=g«8 im, ) 


I 


( COMPUTE ADDRESS IN VERBS WAITISR, 
I ( STALLISR, AND DCORR ) 


590 


4D00 3EA5 

Ml 

c 9 ol 3 E^ 

8 

00 3EBC Cl 
00 3EBD Cl 
00 3EBE Cl 

88 m HI 

AO 3EB1 Cl 
P 3EB2 Cl 
SO 3EB3 Cl 
40 3EB4 Cl 
20 3EBS Cl 
55 3EB6 Cl 

s§ m Hi 

OE 3EB9 Cl 
3B 3EBA Cl 

IIH! HI 

.8 3EC3 Cl 
.4 3EC4 Cl 
2 C! 

C 3EC6 Cl 
08 3EC7 Cl 


1 ( LOHB T 
1 LCIHB 

I 

1 ( PORT A 
( CICT ‘ 

PPC 
DFLAG 
SCT 
DFG 
NLAG 


( ASH ) 


( NLHT ) 


6 « 


8S J88I 8i 


1 


lfADDi8°lSl?I^®|piHIT XSRSET 
TART OF COLLECTION 
CORK 

IRST MINUTE WASTED TO GET MEAN 
i^NUTES e 0 DO 
NT 

I 

I'C« 


STDOS 
! BSA 


ISRRESET 
t 



s'. CR STNRM CR 
* * S. LOOP 


* ' S. . HEX LOOP CR HEX ; 
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l.X MEMOPY USAGE EOR COLLECTION SOPTWARE. 


ADDRESS 

ABBREVIATED 

DEFINITION OF VARIABLE 

IN HEX 

SYMBOL 


80-81 


Tcnp, storage for computation of mult. 

86-87 

RLLL 

Lag table point to LSB 

88-89 

RLL 

Lag table pointer to low byte 

8A-8B 

RLM 

Lag table pointer to med byte 

8C-8D 

RLH 

Lag table pointer to high byte 

8E-8F 

PA 

Contains address for input port 

90-99 


Temp, storage for lag zero values 
during mult. 

A6--A1 

TEMP 

Temporary storage byte 

A2 

HTS 

Height counter used for processing 

A3-A4 

LOLB 

Lag 0 low byte pointer 

A5-A6 

LOAB 

Lag 0 high byte pointer 

A7 

LAG 

Lag counter for processing 

AA 

CICT 

Integration counter 

AB-AC 

LCIHB 

Incoming integration counter high byte 

AD-AE 

LCILB 

Incoming integration counter low byte 

B1-B7 

NLAG 

Offset pointers for computation 

D8-BB 

ASH 

Antenna switching values 

BC 

ANT 

Antenna counter (index to ASH) 

BD 

FPC 

Counter for samples collected vs. EOR 

BE 

DFLAG 

End of record flag (EOR) 

GO 

SCI 

Stall counter for ISR 

Cl 

DFG 

Done flag for sample ready 
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APPENDIX II 

POSTPROCESSING PROGRAMS. 
II. X 3- ANTENNAS. 


REM I CS8A00 

data ii2.7,i3.i;3,4,li2j3il2.3 

i = :? V 

^RINt" " *^‘^(2) ■ .iiCA(3) - .3 

PRINT ; INPUT. PARM FOR OUTPUT PILE* 
PRINT "SLOT#* I ! INPUT SR 
PRINT "VOLUMElSl INPUT VR 
PRINT "INPUT PARM FOR DATA FILES" 

PR INT„" START INC PILE|*>: INPUT MS 
PRINT RENDING FILEI'j: INPUT ME 
PRINT "STARTING VOLUME*>J INPUT VS 
^PRINT^ "ENDING VOLUME''; INPUT VE 
DS ■ CHRS l4) 

PRINT DS;"dPENRESULTS,S*;SR;*,V";VR 
PRINT dSi^WRITERESULTS* 

PRINT ’•ST^T OP FILE* 

ODTM<p n5.»?"T.nQP DPeiiTme* 


INT. D $ ; ^gLggE..RESULTS " 


FOR VZ 
FOR MF 


VE 
MS TO ME 


§ ; -CHRS UfiREM CTRL-D 

fess i 

PRINT DSi»BLOADDRTA'‘|HP»",AS8AOO,SS,V"jVZ 
PRINT " BLOAI)pAT”tHFf',AP9600.S6,V”iVZ 

• ®5fSi§’T0 Si4Sf?t ! 


FOR H - I TO 
Z - 0: FOR N 


sl^iSh2l:S2 


Z * 


10 + PEEK (N) - 176 
35456SS3 » 373768S4 - 38400 
+"256 * (H - 1) ;B2 - S2 + 256 * (H - 1) :B3 - S3 + 


i*OR H- 

B4 » S4 + 128 * (H - 1) 

A - 0; FOR I -.1 TO 3:1A 
A-A+l:L-0: FOR IL - 
L - L + ■ 


(H - 1) :B3 - S3 + 128 * (H - 
O^T(i^io STEP 16 

AC(H7a,LT - 16777216 * PEEK (B1 + IL + lA) + 65536 * PEEK (B2 + IL + 

SliT^e^NEXT^ff N^§^ S 

GOSUB 3000 

|s?rS : • i» - n * 

4 - S4 + 128 - 1) + J3 


B ___ 

FOR IL - 0 TO 80 STEP 16 

¥5’=tRHR-(Smif : fill Hi : sa * 

NfeXT'Ui NEX:^^ft: foil H - 1 TO 8:B1 - SI + 
+ 256 * (H„- 1), + J4:B3 - S3 + 256 * (H - 


1 ) 

STEP 16 


NEXT'iLi next H: foR H - 1 TO 8:B1 
r te.TQjg jTi^ i6:^ - ^ - i 


65536 * PEEK (B2 + IL) + 25 
Jl) * SM(H,J2)) / SQR (AC( 


+ 

B4 

FOR 
CC(H, 
6 * 
CC(H 


256 
S4 + 
IL 
C,L 


256 

1 ) 


(H - 1) + 
J4;L - 6 


J4:B2 - S2 


:C(H,C,L) 


ll?7 

K (B3 + IL) + ,PEE 

i 


EK (B4 + 
122 + L) 


IL) + 6 
- SM(H, 


65536 
Jl) - 


" PEEK (B2 + 
SM(H,J2)) / 


IL) + 25 
SQR (AC( 


2 TO 6 


NEXT ft: 

RESTORE 

• S9.? S TO^§;. for A - 1 TO 3.VFOR L _ 

Nfix| L.4M"HtA)^i^'lU^(iSi“fc^'LiG (AC(H,A,1) )) :AC(H,A,1) 

t NEXT H. _ 

FOR H - 1 TO 8 
FOR L - 2 TQ 6 

acYh.i.l) 


It NEXT A 


FOR H - 1 TO 8 


AC(H,2,L) + AC(H,3,L) ) / 3.0 


IsMX 


1: FOR L - 2 TO 6 

IF MX < AClH.lpL THEN 650 
AC(H,1,L) 


NEXT L 

W < 


MX 

- - . 2 THEN EV(H) - 1 
NEXT H 

REM ST^T OP CURVE FIT 
FOR H - 1 TO 8: FOR C - 1 TO 3 

FOR L - 2 TO 11 
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730 IP NX > CC(H,C,I.) THEN 755 
NX - CC(H,c;i.f 

(K - 6) >3 THEN 920 

1 : 1,8 8 5.: j : f "68t8,s.l 


+ CC(h,C,E - 2) 
+ CC(H,C,K - 1) 


^FOR iio ■ 1 lo^iSo 

*MBXT LQ ' 

CF - CC1h,C,K - 11 - 2 * CC(H,C,K1 + CC(H 

- cc(H,c;x'+'l) 


950 


+Q+Q+4*Pl*X)) 
PI * X) 

,C ,K 4- 1) 


,a,« 1 8i 5 :t t . CP • X • X 

P”ABS'(K - 6) > 4 THEN 925 


NEX. 
REM 
PRINT 
PRINT 
PRINT 
FOR " 
PM 
IP 


i6*)Vx?^‘.n *X) /lO.HX 


OP CURVE PIT 

APPGNpR^gLpiS*>SR;*,V”|VR 


- * APPl 

T • WRITBREl 

.. 5 . 

- .5« GOSUB 7000IT5 -,TT 
__ EV(H)_< > 6 THEN 1430 

**gosub|? 6S1) 

SQR°(TD^H,lf^S TD(H>1) + TT * TT) 
> O.THEN 1430 

SQR (TD(H,2) * TD(H,2) + TO * TT) 


GOSUB '7000 
IP E 


V? ■ 23?'/ SQR (TDTh,3) * TD(H,3) + 

Gosu^jibi* IP p|h[*< > 8*toIn 14; 

Ml • 1§? /, TotHp! ' 


iP 

V3 - 239 
NT » Al 
NT - NT 
GOSUB 
PRINT H , 


TO * TT) 
1430 




- TD 
ABS 


‘>EV(H)',*>SN(H,1)", 


> OTOEN 
NT (VT. ^ 


m 




lT6(tl2) ) + ABS (TO(H,3))) 
lSM(H,2)*,"jSM(H,3)“,’; INT (100 

p*» INT (10 * 


IP EV H) 

PRINT II 
PRINT INT 
PRINT INT 
PRINT INT 
NEXT H 
PRINT^^ CLOSERESULTS" 

NEXT VZ 

PRINT * DONE PROCESSING* 

bS^ hm 

FOR H ■ 1 TO 8t FOR A - 1 TO 3 
I ■ IM(A) 

B2^-^6S536 * PEEK (B1 + I) + 256 * PEEK (B1 + I + 1) + PEEK 

rqvtl -^1 TO 6 , , 

AC|h^A^L) - AC(H,A,L) / (129 - L) - (B2 * B2) / 16384.00 
- Bl^+®2S6» NEXT Ht RETURN 


M5°*^"suB to 0\L PTS OP AUTOCORR 
M then 7030 


IP PM > 

evim - 6 
RETURN 
EV(H) - 5:, return. 
AI « ABS (1 - PM 
K - 1 

FOR L - 2 TO 
IP AI 

Al - 
K - L 
NEXT L 
IP K < 


IF K > 5 THEN 7025 
> 1 THEN llTo 


* T5 


(B1 + 


TT - SQR 
^RETOgN, 

C ■ aIih 


^RBTO 
GH - 
TT ■ ^ 
RETURN 
“EM . 



.4 

a.K) 


(AI / (1 - AC(H,1»2))) * 



-¥+Ntr^cn«(Pt°B,) * .4 


70 


+ AC(H,1,K + D) 


* B - A • 
K - 1 - (p 


i/Jf :ir:1T(ir! : 

Dj' * ( COS T2 * Al) - COS 

— >/ /^. ^ i V.... ^2 * ^1) , 


ii. 


ATN 

; 


^TN T^- 

SSi 

_ COS 
1^^ : K1°S 

/ p3 /^ivi 
AS 7 (1^7 ]v5 
__ TN < 0 THEN TN ■ 

IF BS < AS THEN 7440 
K1 - AStAS - BSsBS - KltTN - TN 
IP AS <_O.OR BS < 0 THEN 7460 


A3) ) - D2 • ( COS 
SIN (2 * A3)) - D2 * 


(2 * Al) 
( SIN (2 


COS 

Al 



1:181 

2 - TN) 

2 - TN) 

iH Him ' »5ii 


+ 1.57079 



: wi * vxji 

: S{ ♦ 





: : !8 


(WY 
TRUE VEL 


XX 


VI 

V2 

V3 


118 

SIN 


SQR (AS) SBX - SQR (BS) 

RN 

» 4 1 RETURN 

1 f 

Ui m\ { M 

* VY{1) + ’W 

„3* Vxiti + W 

sw *^xx - wx * wx 

* XY - WY * WX) / YI :YI 

RETURN 

REM SUB TO CAt APPARENT & 

REH 

'^il4B'7«o “swi'vvlii 

X - - YI / (SP + 1 / SP) 

Y » YI 7 3P'* SP + 1) 

VA - s6r 7X * X + Y » Y) 

GOSUB 9006 ' 

TA ■ TX - TNs IP TA - 0 THEN TA - 

ihi 

T1 >! TXtX - - X:Y ■ - Y; GOSUB 9000 

IF ABS (T1 - TA) < 1.570796 THEN TT - T1 
loTO^iso’^^ - TA) < 1.570796 THEM TT - TX 

Xi..§OSOB 9000:T1 - TX:X - ABS (X) »Y 
GOSUB 9000 

;p ABS (T1 - TA) < 1.570796 THEN TT - T1 
^IP ABS^Itx. TA)..< ;.5.7Q79,6„THEN_TT - TX 

Y 

VC 

RV * SQR ( ABS {(X * X + Y 

VT - VC / SQR (X * X + Y * 

TT ■ TT + TNiTA = TA + TN 

RETURN 

REM SUB FOR ARCTAM 0-360 
TX ■ 0 

IP X - 0.00 THEN 

O.gO THEN ^ 

0 fl^D Y > Q THEN TX 

0 AND Y < 0 THEN TX 

0 AND Y > g THEN TX 

0 AND Y < 0 THEN TX 


1 / (V3 * V3) 




’*li! 


VX 


WX * XY) / YI 


.017 
(TA) ) 


AS) 

THEM 


7990 


- Y 


:P ABS T1 - TA < 1.570796 THEN TT 


Y):^C^-'^*^SQR^1^C) 


IP Y ■ 
IP X > 
IP X > 
IF X < 
IF X < 
RETURN 
IP Y > 
IP Y < 
RETURN 
IP X > 
IF X < 
RETURN 




6^58 

i:ini 


- ATM 
+ ATN 




0.00 THEM 
0.00 THEN 

0.00 THEN 
0.00 THEM 


TX 

TX 

TX 

TX 


1.5707963 

4.71238 

0.00 

3.141592654 
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4- ANTENNAS 


i UK I CSaAOO . 

: SSSi: a ii : J:® 

2 5.2;9^^Vl) - OAa)!OA{12) ■ 2.356194 


jfi ItMT D?|'1ILOADnMA''!Mf^”Aj|8A00''i F 
lil H ZJ TO 8i.?9* i • I toTeV(H. 


a “I»r- 

z ■ Os WI N ■ 35448 TO 35455:2 - Z * 10 ♦ FZK (H) - 176 
111 NEXT R 

290 81 - 35328:82 - 35456:83 - 37376:84 - 38400 
210 FOE H " 1 TO 8 
215 11 - 81 + 256 * (H - 1):B2 - 82 ♦ 256 * <H - 1):B3 - 83 + 128 * (H - 

2^9 14-84 ♦ 128 * (H - 1) 

5: FOE U - 0 TO 15 STEP 5 


216 14 -84 * 1 
220 A - 5: FOE 
225 A - A - 1 :L 


128 * (H - 1) 

U - 0 TO 15 STEP 5 
L - 0: FOE IL - 0 TO 80 STEP 16 


230 L - L 4- 1 

235 4C(H,A,l 5 - 16777216 * PEEK (B1 ♦ IL ♦ lA) 4- 6 5 536 * PEEK (B2 ♦ IL 

l« Slir'Ai.r'S? SSrV’ * “> * 

244 GOSDB 3000 

245 FOE C - 1 TO 6; READ J1,J2.J3,J4 

250 FOE H - 1 TO 8:81-81 * lib * (H - 1) ♦ J3;B2 - 82 ♦ 256 * (H - 1) 
... S3 + 128 * (H - 1) f J3:L - 5 


SM(H,J2)) / 8QE (AC( 


. J3:B3 - S3 + 128 * (H - 1) f J3:L - 5 

251 B4 - 84 4' 128 * (H - 1) * J3 
255 FOE IL - 0 TO 80 STEP 16 
260 L — L ♦ 1 

265 CC(H,C^) - 16777216 * PEEK (B1 + IL) + 65536 * PEEK (B2 + IL) 25 
6 * PEEK (B3 + IL) f PEEK (B4 + IL) 

270 / <134 - L) - 8M(H,J1) * 8M(H,J2)) / SQR (AC< 

275 Nixr’ll.: NECT 4: #OE,H - 1 TO 8:B1 - 81 ♦ 256 * (H - 1) ♦ J4:B2 - 82 

J7. r»vi! : ,v{i« :*» * • !h - 15 • «.L . B 

285 CpH.CjL) i^l6777216^ PEEK (B1 1 IL) + 65536 • PEEK (B2 ♦ IL) ♦ 25 
_ 6 * tEEK (B3 + IL) t PEEK (B4 ♦ IL) 

295 ! ’122 + L) - SM(H,J1) * SM(H,J2)) / 8QE (AC< 

296 HiXT’lL: NEXT 4: 4 eKT C 
300 RESTORE 

To 8: rOE A - 1 TO 4: FOE L - 2 TO 6 
AC(H,A,L) - AC(H,A,L) / AC(H,A,1) 

NEXT i:AC(H.A.l) ~ 1: HEXI A;'reXT H 
FOR H - 1 TO 8 
FOE L - 2 TO 6 

, 0 SI - AC(H,1,L) + AC(H,2,L) + AC(H,3,L) 

510 82 - ACCH.I L) ♦ AC<H,2 L) + AC(H 4.L) 

520 S3 - AC(H,2,L) + AC{H,3,L) ♦ AC(h;4,L) 

>acU:4.lJ 

ii8 18 : ll ^ ^.8 

590 FOR U - 1 TO 8: FOR 8 - 1 TO 4 
600 K - 1:MX - 1: FOE L - 2 TO 6 
610 IF MX < AC<H,S,L) THEN 650 
620 MX - AC(H,S,l) 

630 K - L ’ ’ 

640 NEXT L 
650 AH(H.S) - MX 
660 IF K < 3 THEN EV(H,S) - 1 
NEXT 8: NEXT H ’ 

Em START OF CDEVE FIT 
CA(1) - .1:CA(2) - .2:a(3) - .3 
696 CA<4) - ,1:CA(5) - .2:a(6) - 1 
700 FOR H- 1 TO 8: FOE C - 1 TO 6 
795 IF EV(H,1) < > 0 THEN 920 

s “fo.- f'iVtj’if : ‘ 

;i8 

750 K - L ’ ’ 

755 NEXT L 

760 IF ABS (K - 6) > 3 THEN 920 

■" 85 : 88ig;8:5 : 1 : 

CH - CC(H,C,K + 1) - CC(h|c’k - 1> 

ss s : sKvci 7 1: ■ ■’ 

5 : : f.-.'/*.* ' *•= 


THEN 755 


* CC(H,C,K > 2) 



72 


fP! 


i - S > X ♦» + X*(Q+Q*Q + 

IF MS (Fi)_<:^.090i thw'mo ^ ^ [ 
2 t * 9 ♦ 12 * PI * X) 


♦ 4 * M • X)) 


m x^- X'- Fx / ft 

mO next.uj 

920 P(H,C) ■ - 2:TD(H.C) - 0: GOTO 9S0 

a Mi-.'ifl'i ; d<«‘ ♦ « MO . 

fiS SSi‘'"en» 

RSI BID or CURVE FXT 
?R H ?. 1 TO ■ 

TD' 

TO 
TO 
TO 
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l|g Wi*’ ■ 

[I|0 PRINT "TIME-"iZ 
i200- FOR H - 1 TO ft 
2lg FOR S,- 1 TO 4 



l220 IF EV(H.S) 

PM - FT''^ - 

G08UB 


?P EV(H,sy < 


< > 0 THEM 1430 


RNTll?Hfs,it^2 TQ(H,8,1) ♦ TI • TI) 


^EN X430 , 

H,S,2) * T(}(H,8,2) + TI * TT) 


9.S) < >,0 THEM I46S 

fe HiMrSi fSjMi'iSiS’Mj) 

SlI ^7.246'* TM)\«; !«' 


TO 4 
+ 256 


i260 VI - X69 7 SO 
1270 PM - PQ(H,S,2; 

1280 COSUB 70flQ 
X290 IF EV(H,S) < > 0 

1300 V2 - 16?/ SQR 
1310 PM - PQ(H.S,3J 
1320 GOSUB 7O0O 

ilfo SVPlts «l(H,S,3) ♦ TT * Tt) 

Ilfs t 

ins '™ ^ .5: GOSUB 7000:T5 - TT 
1390 Vl°S”?65^?°i«R ^ ° 

iWlM : ill 'i iliiif, 

Hi! ^ 

1430 PRINT H‘,S'/,"5EV(H,8) 

1431 IF EV(hIs) " " 

1432 PRINT 

1433 PRINT 

1434 PRINT v.w/ , . 

1435 PRINT IHT (10 * RV) 

1465 NEXT S: NEXT H 
1520 NEXT MF 
2000 NEW : END 
3000 B1 - 35572 

3010 FOR,U - 1 TO 8: FOR A - 1 
3020 I - (A - 1) * 3 
3030 B2 - 65536 * PEEK (B1 * I) 

1 + 2} 

3040 FOR L -7 TO 6 

3060 mxtV^ *' 

3065 IF ACIH.A.l) ,< 0 THEN GOSUB 4000 

3070 SM(H,A} - fc2 / 128 

3080 NEXT A;B1 • BI + 256: NEXT H: RETURN 

8!8 SaVS.) I'W* ■ ’ 

4020 RETORN 

6999 RBI SUB TO CAL PTS OF AUTOCORR 

7000 IF PM > - AMfH.S) THEN 7030 
— 1 EV(H.S) - 10 

) RETORN 

.,.3 EV<H,S) - 11: RETURN 
7030 AI ABS (1 - PM) 

7040 K - 1 

7050 FOR L - 2 TO 6 

7060 IF AI < ^ (PM - AC(H,S,L)) THEN 7090 
7070 AI - ABS (PM - AC(H,S,l)) 

7080 K * L 

7090 NEXT L: IF K > 5 THEM 7025 

Ills • •* 

;{|§ J :■ 

7160 IF MS7Ay> .0001 THEN 7190 

7170 TT - (K - 1 ♦ (PM - C) / (2 * B)) * .4 

7180 RETURN 

7190 GH - B * B - A * (C - PM): IF CH < 0.0 THEN 7170 


AB8 (TQ(H,S,3))) 


do’* NT)","; 


PEEK (B1 + I + 1) + PEEK (B1 + 
(B2 * B2) / 16384.00 


a 




,MH , SU|„TO aL.StUKE,nUUI , 

sH 1, ? : : 8 : ? 

I? £ ?i^£ ' « • **) - COM* 


,5 * a ♦ cos d * ai - 
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* Ai)) - D 2 * < col (J * Al) - COI 
ilH (2 • AS)) - D 2 • ( IXN (2 * Al 


* (I ♦ COS 12 * 

* (1 - COS (2 * 

* K2 - ri « ik 


Al - 1 M 
Al > 2M 
A2 - 1M 
A2 - w; 


•s : p 1 s\i> - II 

™Vw*+ jiUi 
ir BS < AS IMn 74A0 


ir BS < AS 


„ ,V1 # VI) - Kl „ ... 

I <v2 * n) - tA I In * nj 

Ml TO - TO ♦ JiUlM 
Hn 7440 

»§'•* ; KiiPf.r.P* ♦ 1.5707J 


(V2 * T2)) 
Jvi * VI)) 


Kl - ASiAS - BStlS - KliTO - T 
!r AS < O.OK.BS < 0 TODI,74(0 


RITVBM 

4»[S,S5 


-^9« BETUim 


, SLplE riT FOR APFUITOT VEL 

s : m J t 

* mi) ♦ W2 * VX(2) 


- I / (V 3 * V 3 ) 


W1 * VX(1) * VX(1) 
Wl ♦ W2 ♦ W3 

f n * XX - wx * wx 
SV * XT - WT • W 


* VX{3) * VY(3) 

• VX(3) • VX(3) 


RETUTO * ^ “ W • *•*) / ■ <Wr ^ M - W * *») / U 

RTO SUB TO CAL APFARBIT 4 TKII VEX. 

X(l) • VI * COS (A1)!VY(1) - VI * SXM (Al) 

X(2) - VJ * COS (A2)sVY(2) ** V2 • SXH (A2) 

?i3L“ YL* COS (A3)IVYU) - V3 * SXH (A3) 


VX(1) • VI * COS (A1)!VY(1) - VI * SXM (Al) 

VX(2) - V2 * COS (A2)sVY(2) ** V2 • SXH (A2) 

'coiuB’7480 (a3)iVY(3) - V3 * SXM (a3) 

X ■ - YI"7 (8P ♦ 1 / 8P) 

Y ■ YX / (fep * 8P ♦ 1) ' 

VA - SOR (X * X + Y * Y) 

GOSUB MOO 

TA - TX - TOi IF TA- 0 TOIM U ■ .01? 

n ■ VA * ^ SI^(TAt - 3P * 008 (TA)) 

X ■ AB^ (S?lf * A8^/*^SQR*fBS ^|p * SP * AS) 

“!'*»» / (Y • AS)) < .0001 THEM 79« 

GOSUB 9000 

T1 - TXsX ■ - XiY ■ - Yi GOSUB 9000 
IF ABS (T1 - TA) < 1.570796 TOEMTT - Tl 
^ 1.5707 96 TOEM TT - IX 

GOTO 8030 

K:^r.„Z 54„S°SUB 9000 !T1 - TX;X - A3S (X>!Y - 
GOSUB 9000 

IP ABS (Tl - TA) < 1.570796 THEM TT - Tl 
RS_(tX - TA) < 1.570796 TOW TT - TX 


9000 !T1 - TX;X - ABS (X>!Y - 

GOSUB 9000 

ir ABS (Tl - TA) < 1.570796 TOW TT - Tl 

IF ABS (tX - TA) < 1.570796 TOW TT - TX 

XO - XiYO - Y!X - YI / (YO / XO - 3P) 

Y - (YI * YO / XO) / (TO / XO - SP) 

VC - XO * XO V YO * YO 
RV - SQR ((X * X ♦ Y * Y) / VC - 1) 

VT - V^/ SQR (X * X ♦ Y * Y):VC - SQR ( 

TT ■ TT ♦ TWiTA - TA + TO ^ 

RETURN 

REH SUB FOR ARGTAN 0-360 
TX ■ 0 

IF X - 0.00 TOW 9060 
IF Y - 0.00 TOW 9090 
iE 5 ? S ^ > 0 THD* T* " ATN (Y / X) 

IF X > 0 AND Y < 0 THEN TX - 6.28 - ATO 

IF X < 0 AND Y > 0 TOW TX - 3.1415 - AT 

IF X < 0 AND Y < 0 TOW TX - 3.1415 ♦ AT 

RETURN 

IF Y > 0.00 TOW TX - 1.5707963 
IF Y < 0.00 TOW TX - 4.71238 
RETURN 

IF X > 0.00 TOW TX - 0.00 

TX - 3.1415926 54 

RETURN 


SQR (VC) 


SA < ! < f 

♦ ATO (y / X) 
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APPENDIX III 
CVBER PLOTTING PROGRAMS. 


nociAK sc(iMm,ouTfuT,cm,Tm3H:sno 
nOGRAM 8C IS USED TO nODOCE SCATTBt 
riots THAT COMPAIE C.S. WDIDS DATA TO 
P.R. UIMD DATA. SEQUUn A PIU CSFR 
TEAT CONTAINS TNE WIND VALUU FOR 
" NAXI " VALUES AT EACI ALTITUDE. WILL 
ASE rOR THE DATE | tWNni,0AY 
EnEINALS: EETA / GUB.ZETA 
PLOT COMMAND! PLOTZ.TAPeiS/LEMCTE-IO/TIME-ZA 
WEITTIN BY RAYMOND RUGGERIO 19R2 
DIWNSION CS(7 ,20) ,PR(7 ,20) ,R(20) , Y(20) ,RP( S) 

DATA »r/70.5,75.,79.5,S4.,eE.>,93.,»7.5,I02./ 

PRINT 90 

90 POBMAT(''INPUTiHAXI,TMTt,TDAY") 

READ *,HAn,TMTH,TDAY 

RRAD(3,40)((CS(I,J),PR(I,J),J-1,MAXI), 1-1,7) 

40 PORMAT(2P10.1) 

CALL PL0T8(0.,0.,26) 
call PL0T(l.0,3.0,-J) 

DO 150 IH-l,7 
DO 120 M-1,HAX1 
X(M)-PR(IH,M) 

Y{M)-CS(IH,M) 

120 CONTINUE 

51- 0 

52- 0 

551- O 
SXl-0 

DO 950 H2-1,NA» 

51- SUFR(IH,H2) 

52- S24CS(IH,M2) 

8SI-SS1+PR(IH,M2)*PR(IH,M2) 

552- S82^S(IH,H2)*CS(IH,H2) 

950 SX1-9X1+CS(IH,M2)*PR(IH,K2) 

31-81/ FLOAT(HAXI) 

S2-82/FLOAT(HAXI) 

8S1-SS1/FL0AT(MAXI) 

S82-8S2/FLOAT(MAU) 

SXI-8X1/FL0AT(MAXI) 

FHO-(SXl-Sl*S2)/8qRT((S8l-8l*Sl)*(SS2-S2*82)) 

PRINT *,IH,PH0 

CALL AXIS(0.,0.,"PR VELOCITY",-!! ,4. ,0. ,-100. , 50.) 
aiL AXIS(0.,0.,"CS VEL0CITY",!!,4.,90.,-!00.,50.) 
aiL 3YMBOL(2.0,4.5!,.!05,"ALT- ”,0.,5) 
aiL m!MBER(2.42,4.5l,.!05,HP(IH),0.,!) 
aiL 8WBOL(2.84,4.51,.!05," XM",0.,3) 
call 8WB0L(2.,4.72,.!05 ."DATE" ,0 . , 5 ) 
aiL miMBER(3.1,4.72,.!05,lMra,0.,-!) 
aiL SWB0L(3,3!,4.72,.!05,"/",0.,!) 

CALL WJMBER<3.4!5,4.72,.!05,TDAY,0.,-!) 
aiL SYMbOL<3.625,4.72,.!05,"/82",0.,3) 
call 8YMB0L(2.,4.3,.105,"PH0-",0.,4) 
call NUHBER(2.42,4.3,.!05,PHO,0.,3) 

CALL PLOT(0.,0.,3) 
call PL0T(4.,4.,2) 
aiL PLOT(4.,4.,3) 
au PL0T(2.,2.,-3) 

DO 130 >!,HAXI 

!30 CALL SWBOL(X(J)/50.-.054,Y(J)/50.-.054,.U,'V',0.,!) 
!50 aLL PLOT(6.5,-2.,-3) 

CALL PLOT(0.,0.,999) 

STOP 

END 
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nOGMM REC0V(INm,0UTm,KIN,10UT,TArB3-t0DT,TAm->XN, 
♦RCS,TAPE4-RqS) 

noCRAM RECOV IS USED TO RROCEBS DATU 
TRAHSEERKED TO THE CYBER EROH THE ATPLE COHPUTER 
REQUIRES INPUT PILE TO HAVE THE 
NAM RIN WHICH CAN BE ACCONPLIBED BY 
USING THE RENAME CMOS RENAME, RIN-PI LEE 
PRODUCES OUTPUT PILES ROOT AMD RCS 
ROUT IS THE INPUT DATA REPORMATTED AND 8CREOIED 
RCS IS THE MEAN UIMO VALUES PROJECTED 
ONTO C.S. LINE OP SIGHT 
TAPE99 IS THE PLOT COMMAND PIU PRODUCED 
IT CONTAINS THE IMPOEMATION TO PRODUCE THE 
PLOTS POR NS, EH HIND AMD CORREUTXON TIMl AND POHU 
EXTERiVALSs ZETA / GRAS,ZETA 

PLOT PILE USING CMDs PLOTZ,TAPB99/LENG1H-170/TIME-30 


WRITTEN BY RAY RUGGERIO I98Z 


DIMENSION VT( 8,2AO) ,TT( 6.240) ,AX( 8,240) ,BX( 8,240) , 
fIN( 8,240 ) ,EW( 8,240) ,AR( 8,240) ,VPV(8) ,RNS( 8,240 ) 

AH( 8,2), HD(S), PA1(8,240),PA2( 8,240), PA3( 8, 240), TS( 8,240) 
+,VCS( 8,240), VCSM(O) 

INTEGER. TINE 

COMMON VT,TT,AX,BX,TN,RMS,EH,AR,AH 
aLL PLOT5(0.,0.,99) 
aLL PLOTa.,1.,-3) 

RVT«200 
RNTT-3 
VMAX-200 
T5MIN-15.0 
DO 50 HI-1,240 

READd , 10 ) (HD( I ) , I-l ,4) ,TIME 
IP(EOPa),HE.O)CO TO 900 
10 P0RHAT(4A1,1X,I7) 

IP(HD(1).HE."T")C0 TO 800 

IP(MI.Eq.l)aLL trcv(time,tmth,tday,thr,tmin) 

DO SO IH-1,8 

READd ,*)M, IE ,PA1(H,HI) ,PA2<M,MI) ,PA3(N,HI) .T5(H,MI) 
IP(IE.NE.O)GO to 60 

REAO( 1 ,*)VT(H.HI) ,TT(M,MI) ,AX(M,HI) ,BX(N,HI) ,VC, 
tTM(M,MI),lNT,RV 
IP(RV.GT.RVT)GO TO 60 
1P(RNT.GT.RMTT)G0 TO 60 

ip(vt(n,mi).gt.vhax)go to 60 
IP(T5(M,MI).LT.T5HTN)00 TO 60 
GO TO SO 
60 VT(M,MI)-0 
TT(M,MI)-0 
AX(M,MI)-0 
BX(M,MI)-0 
TN(M,MI)-0 

SO IP(T5(H,MI).LT,0.0)T5(M,MI)-0.0 
900 MAXX-HI-1 

1P(HAXI.LT.2)C0 to 800 
DO 200 M-1,HAX1 
DO 200 IH-1,8 

RNS( IH ,H) — VT( IH ,H)*COS( .0174533*TT(IH ,M) ) 

EW( IB ,M)— VT(IH ,H)*SIN( .0174533*TT( IH ,M) ) 
AR(IH,M)-BX(IH,N)/10.0 
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200 CONTIMUE 
DO SOO 

ll>0 

Cl-0 

DO 920 K-l.)Un 
VCI(»,M)-UI(»,M)*.9I77*-W(»,M)*.009 
ir(TT(»,M).B().0.)CO 10 520 

ci-cui 

920 COMTINDE 

VCIM(IH)-SI/C1 

500 COMTINDE 
nitE(4,540)TIME 

540 r0i»UT(lHl,"Tn«-'*,I7) 

Wim(4,550)«VCI(II,N),N«i,Nai},IH-l,8> 

590 70mAT(6F20.1) 

WEm(4,555)<VC8M(J),J-l,8) 

555 rOIIUK//" MEAN VA10E87(r20.1)) 

aiL ADISP(T5,KAXI,1NTU,TDAT,TU,1MIM) 

CALL mSF(PAl,rA2.PA3,IIAZl,1Nn,TDAT,TU,'»lM) 
Vim(3,450) 

450 rOEMAT(lHl,2:;;”HIK,VT,Tr,HS,EH,a,TM,R>Nn.lULr TZME**) 
DO 300 IH-I,8 
DO 310 M-l,HAn 

WEITE(3,410)M,VT(II,M).TT(X8,M),BM8<XI,M)^EH(XI,H) 

4^ ,At( lU ,H) iTNaa .M) , PAK XM ,K) ,15(11 .M) 

4X0 P0BMAT(2X,X6,8P15.4) 

310 CONTINUE 
WRXTE(3,420) 

420 FOEMAiaHl) 

300 CONTINUE 

CAU AVER(MAXX) 

PRINT SOI 

501 F0EMAT( "INPUT VECTOR FOR DESIRED PLOTS") 

PRINT $02 

502 FOINAT(" A 1 MEANS PLOT, A 0 MEANS DO NOT PLOT") 

READ *,(VPV(J),J-1,8) 

DO 987 ID-1,8 

IF(VPV(IU).Eq.l.O)aLL DlSPU(KAn,IO,TMTB,TDAT,nR,1MlN) 
987 CONTINUE 

CALL PLOT(0.,0.0,999) 

CO TO 911 
800 vniTE(3,430) 

430 FORMATC SEQUENCE ERROR") 

911 STOP 
END 

SUBROUTINE ADISP(T5 ,MAn,1NTH,TDAY,1U,TNIN) 

DIMENSION T5(8.240) 

DO 10 J-1,8 

DO 10 M-1,MAXI 

T5(J,M)-T5(J,M)/100. 

10 IF(T5(j,M).GT.2.0)T5<J,M)-0.0 

CALL AXIS(0.,0.,15HTIME IN MINUTES,-15,6.,0.|0.,10.) 
aLL AXIS(0.,0.,13HALT1TUDE (EM), 13, 7. ,90. ,70. 5, 4.5) 
ESC- 10. 

Y8C-1.5 

aLL STMB0L(3., 8.3, .105, "START TIME-", 0.0, 11) 
aLL NUMB£R(4.155,8.3,.105,THR,0.,-1) 
aLL SYMSOL(4.365 .8.3 ,.105,":", 0.0,1) 
aLL HUMBER(4.47,8.3,.105,THIM,0.0,-1) 
aLL 8YMBOL(3., 8,51, .109, "DATE ”,0.0,5) 
aLL NUMBER(4.1,8.51..105,1}ITH,0.0,>1) 
aLL SYMBOL(4.31,8.$1,.105,"/",0.0,1) 
aLL NUMBER(4.415, 8.51,. 105 ,TaAY,0.0,-’l) 
aLL SIMB0L(4.65, 8.51, .105,”/8Z” ,0.0,3) 
aLL PLOT(-.0625,8.,3) 
aLL PL0T(.0625,8.,2) 
aLL PL0T(0.,8.,2) 
aLL FL0T(0.,9.,Z) 
aLL PLOT(-.0625,9.,2) 
aLL PL0T(.0625,9.,2) 
aLL PL0T(.0625,9.,3) 

aLL SWB0L(. 25, 8. 51,. 105, "CORRELATION TIME",0.,16) 
aLL SWBOL(.25,8.3,.105,"1.5 SEC”,0.,7) 

DO 200 lR-1,8 


V 
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CALL PL0T(6.,0.,2) 

CALL rLOT(6.,0.,3) 

CALL r/^Ta/XSC,T5(IK,l)/T8C,3) 
m 21( v)'>2,MAXI 

2i0 ir(T5v;gvJ}«KK.0.)CALL PL0T(rL0AT(j)/XSC,TS(XI,J)/TIC,2) 
aLL PL0T(rLDAt(HAn)/ZSC,TS(n.HAn)/YIC,3) 

CALL PL0T(0..1.,-3) 

200 CONTINUE 

aLL PWT(8.5,-8.,-3) 

RBTUIN 

END 

SUBkOOTINE FDISP(FA1,PA2.PA3,HAXI,TIITH. that, TRK, WIN) 

DIMENSION PAU 8, 240 ) , PA2( 8, 240 ) , PA3( 8, 240 ) 

DO 10 IH-1,8 
DO 20 J-l.MAXI 

PA1(IH,J)-10*ALOG10(EXP(PAU1H,J)/100.)) 

IP(PAl(IH,J).EQ.0.)G0TO 40 
CO TO 20 

40 ir(FA2(XH,J).NE.0)PAl(lH,J)-10*ALOGl0(EEP( 

+PA2(IH,J)/100,)) 

IP(PA1(IH,J).EQ.O.)PA1(IH,J)-10*ALOG10(EXP( 

•t-PA3(IH,J)/l00.)) 

IP(PA1(IH,J).EQ.O.}PA1(IH,J)-SO.O 
20 CONTINUE 
10 CONTINUE 

8PMIN-PA1(1,1) 

SPMAX-PAl(l.l) 

DO SO IH-1,8 
DO 50 >1,HAXI 

IF( PAl < IH , J ) . GT. 3PMAX) SPMAX-PAl ( IH , J ) 
so IP(SPMIN.GT.PA1(IH,J))SPMIN-PA1(IH,J) 

DO 60 IH-1,8 
DO 60 J-1,HAXI 
PAl < IH , J ) -PAl ( IH , J ) -SPMIN 
60 CONTINUE 

aLL AXIS(0.,0.,15HTIME in MINOTES,-15,6.,0.,0,,10.) 
aLL AXIS(0.,0.,13HALTITUDI (iCM), 13, 7. ,90. ,70. 5, 4.5) 

XSC-IO. 

YSC-FLOAT(IFIX( SPMAX) ) 

aLL STMB0L(3..8.3..105."START TIME-". 0. ,11) 
aLL NUMBER(4.155,8.3,.10S,THK,0.,-1) 
aLL SWB0L(4.365,8.3,.105,":",0.0,1) 
aLL NUMBER(4.47,8.3,.10S,'miN,0.,-l) 

S.ALL SWB0L(3.,8.51,.105,"DATE ", 0.0,5) 

aLL NUMBER(4.1,8.51,.105,TMTH,0.0,-1) 

aLL STMB0L(4.31,8.51,.105,"/",0.0,1) 

aLL NUMBER(4.415,8.51,.105,TDAY,0.,-1) 

aLL SYMBOL(4.65,8.51,.105,"/82",0.0,3) 

aLL PL0T(-.0625,8.,3) 

aLL PLOT<.0625,8.,2) 

aLL PL0T(0.,8.,2) 

aLL PLOT(0.,9.,2) 

aLL PLOT(-.0625,9.,2) 

aLL PLOT(.0626,9.,2) 

aLL PLOT(.0626,9.,3) 

aLL SYMBOL( .25, 8.51, .105, "POHEB" ,0.0,5) 

aLL NUMBER(.25,8.3,.105,YSC,0.,-1) 

aLL SYMBOL(.25,8.09,.105,*’OB",0.,2) 

DO 200 IH-1,8 
aLL PLOT(0.,0.,3) 
aLL FLOT(6.,0.,2) 
aLL PL0T(6.,0.,3) 

CALL PLOTa./XSC,PAl(IH,l)/Y8C,3) 

DO 210 J-2,HAXI ' 

210 aLL FL0T(FL0AT(J)/XSC,PA1(IH.J)/YSC,2) 

aLL PL0T(FL0AT(MAZI)/X8C,PA1(IH,HAXI)/Y3C,3) 
aLL PLOT(0.,1.,-3) 

200 CONTINUE 

aLL PL0T(8.5,-8.0,-3) 

RETURN 

END 

SUBROUTINE WCTC TIME, WTH,TDAT, THE, WIN) 

INTEGER TIME 

TMTH-IFIX(FLOAT(TIME ) / lOOOOOO . ) 


inAT-mx(rLOAT(Tim)/ioooo.)-'DfTH*ioo. 78 

m-mx(n.oAT(TmE )/ 100 . )-mrH*ioooo ,-mY*ioo . 

mXX-n.0AT(Tna ) - 100*7M- 10000*mY-1000000*TMTH 

UTUM 

END 

SUIROUTINE DISPW2(X,Y.MAXI,WD,I,Timr,TDAY,THX,TMIN) 

DIHBM8I0H X(240),Y(240),HP(8) 

DATA HP/70. 5, 75., 79.5,84., 88.5, 93. ,97. 5, 102./ 

ORIGINAL PAQE IS CAJ-I* AXIS(0.,0.,15HTim W MIKOTE8,-15,6.,0.,0.,10.) 

_ _ rvilAI itV ip(wd.bq.o.)call axis(o.,o., 

OF POOR QUAUnf t'lWXTH TO SOUTH VELOCITY (M/8)",29,8.,90.,-200.,50,) 

IP(WD.EQ.l.)aLL AXIS(0.,0.,"EAST TO HE8T VELOCITY (M/8)", 

♦27, 8.0, 90, ,-200. ,50.) 

X8C-10. 

Y8C-50. 

aLL SYM»«L(3.,8.3,.105,"ALT- ",0.0,5) 

MLL NUMBER(3.42,S.3,.105,UP(I),0.,1) 

CALL SYMB0L(3. 84,8.3, .105," KM",b.,3) 

aiL SYMBOL(3.,8.5l,.105,"STAXT TIME-", 0.0, 11) 

CALL NUMBER(4.155,8.51,.105,THR,0.,-1) 

CALL SYMB0L(4.365, 8.51, .105,":", 0.0,1) 

CALL NUMBEH(4.47, 8.51,. 105, TMIN,0. 0,-1) 

CALL SYMBOL(3.0,8.72,.105,"DATE ",0.0,5) 

CALL NUMBEH(4.1.8.72..105.TMTH,0.0,-1) 

CALL SYMBOL(4.31,8.72,.105,"/",0.0,1) 

CALL MUMBER(4.415, 8.72,. 105 ,TDAY,0. 0,-1) 

CALL 3YMB0L(4.65, 8.72, .105, ’782", 0.0,3) 
aLL PL0T(0. 0,4. 0,-3) 

DO 10 J-1,HAXI 

ir(r(j).BQ. 0 . 0)00 to lo 

aii SYMB0L(X(J)/X8C-.054,Y(J)/Y8C-.054,. 11, "+",0.0,1) 

10 CONTINUE 

aLL PL0T(8.5,-4.0,-3) 

RETURN 

END 

SUBROUTINE AVER(MAXI) 

DIMENSION VT(8,240),TT(8,240),AX( 8,240), TM(8, 240), 

+RNS( 8,240 ) ,EW( 8,240) ,AR( 8,240) ,AH( 8,2) ,BX( 8,240) 

COMMON VT,TT,AX,BX,IM,RMS,EH,AR,AH 
INTEGER H,SFF(8'I 
DO 10 H-1,8 
81-0 
82-0 
Cl-0.00 

DO 20 M-1,MAXI 
IF(VT(H,M).BQ.O)CO TO 20 
Cl-Cl+1 

51- 81+RNS(H,M) 

52- S2+EW(H,M) 

20 CONTINUE 

SPF(H)-C1 
AW(H,1)-S1/C1 
10 AH(H,2)-S2/C1 

HRITE(3 ,700)(AH(H,1) ,AH(H,2) ,SFF(H) ,H-1 ,8) 

PRINT 700,(AH(H,1) ,AH(H,2) ,SPF(H) ,B-1,8) 

700 F0RMAT(1H1//2X," MEAN WIND VALDES NS AMD EW "/(2F10.1,5X,I5)) 

RETURN 

END 

SUBROUTINE DISPW(MAXI,H,TMTH,TDAY,IHR,TM1N) 

DIMENSION X(240) ,Y(240) ,VT( 8,240) ,TT( 8,240) ,AX( 8,240) , 

+BX( 8,240 ) ,TM( 8,240) ,RNS( 8,240) ,EH( 8,240) ,AR(8, 240) ,AH( 8,2) 

COMMON VT,TT,AX,BX,TH,RMS,EH,AR,AM 
INTEGER H 
DO 10 J-1.MAXI 
X(J)-J 

10 Y(J)-BNS(H,J) 

HD-0.0 

aLL MSPW2(X,Y,MAXl,HD,H,TMTH,TDAY,TaR,TMIH) 

DO 20 J-1,MAXI 
20 Y(J)-EW(H,J) 

HD-1.0 

aLL DISPH2(X,Y,MAXI,HD,H,TMTH,TEAY,THR,TMIN) 

RETURN 

END 
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APPENDIX IV 

OPERATION MANUAL FOR PARTIAL-REFLECTION SYSTEM. 


This is a step-by-step procedure that must be followed to operate the 
part ial-rcf lection-drifts experiment. 

1) Connect the pulser to the radar controller to obtain the 400 He 
control signal. 

2) Turn on the power switches on the pulser, A/D, receiver, relay 
switching network, and all ac. switches on the transmitter panel. 

3) Unground the T/X-antenna and the R/X-antenna systems. 

4) Slowly turn the filament voltage for the final stage amplifier 
tubes to 5.5 volts. 

5) Turn the bias voltage to 28 volts. 

6) Turn the switch located under the high voltage (5 kV) variac on, 
the bias voltage will jump to 35 volts. 

7) Turn the variac clockwise until the high voltage is up to 2000 
volts. 

8) Press the low voltage reset on the 2 kV power supply. 

9) Slowly turn up the 2 kV power supply up to 1100 volts while 
maintaining 2 kV on the high voltage power supply. 

10) Turn the high voltage up to 3000 volts. 

11) The T/X and R/X are now on the air and a'l that is required no;-; 
is to bring up the computer system. 

Computer System: ^ 

1) Place the following peripherals in the appropriate slots of the 


computer; 
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CorvuB disk drive — — — — slot #6 

Parallel interface Slot #7 

Arith proc Slot #5 


2) Sign on the system as user "PR". 

3) Type "BRUNDRIFTC0L.S6 ,V49". 

4) Type "MAINDRIFT". 

5) The program will start putting data on volume 30 and continue to 
volume 67. File numbers will range from 0-19 for each volume. 

6) Collection computer prompts the volume# and file# that corresponds 
to the last minute of data that was collected and stored. 

To bring up the processing computer; 

1) Sign on as user "PR". 

2) Type •'MAXFILES 1". 

3) Type "BRUNBIN.DRIFTL2 ,S6,V49". 

4) Enter information that is requested. 

Example; 

INPUT FARM FOR OUTPUT FILE 
SL0T#4 
VOLUME #254 

INPUT FARM FOR DATA FILES 
STARTING FILE# - 0 
ENDING FILE# =19 
STARTING VOLUME # - 50 
.ENDING VOLUME# = 52 

Program will display what file and volume it is working on. 

3) The text output file generated will have the filename of "RESULTS". 
6) The order of the results printed for each minute is: 
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LINE#1: Date and Time 

*LINE#2; ALT#, Error vector, Power of antenna 1, Power antenna 2, 

Power antenna 3, Autocorrelation halftime times 100. 

LINE#3; True velocity. Direction of drift in degrees. Semimajor 
axis "a". Axial ratio, VC , Rotation angle of ellipse, MTD, Random 
variation times 10. 

Once the data are stored on the floppy disk the results can be trans- 
ferred to the University computer for plotting via a modem. Plotting 
routines for the computer are available for this. 

*If the error vector is not equal to zero the next altitude ie printed. 

The error vector indicates that the drift analysis failed to generate a 
true velocity. 
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